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Abstract 

We study general models of holographic superconductivity parametrized by four arbitrary func- 
tions of a neutral scalar field of the bulk theory. The models can accommodate several features of real 
superconductors, like arbitrary critical temperatures and critical exponents in a certain range, and 
perhaps impurities or boundary or thickness effects. We find analytical expressions for the critical 
exponents of the general model and show that they satisfy the Rushbrooke identity. An important 
subclass of models exhibit second order phase transitions. A study of the specific heat shows that 
general models can also describe holographic superconductors undergoing first, second and third (or 
higher) order phase transitions. We discuss how small deformations of the HHH model can lead to the 
appearance of resonance peaks in the conductivity, which increase in number and become narrower 
as the temperature is gradually decreased, without the need for tuning mass of the scalar to be close 
to the Breitenlohner-Freedman bound. Finally, we investigate the inclusion of a generalized "theta 
term" producing Hall effect without magnetic field. 
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1 Introduction 

The AdS/CFT correspondence [U [21 [3] has been recently applied to provide a gravitational description 
of systems that undergo superfluid or superconducting phase transitions. Holographic superconductors 
[4l[5l[6] reproduce various features of the physics of BCS superconductors, and they represent a promising 
avenue for constructing tractable models of strongly interacting systems like high superconductors or 
other systems exhibiting non-Fermi liquid behavior. 

The basic model of [5j (henceforth, the HHH model) depends only on one free parameter and obviously 
does not have enough room to incorporate non- universal features of the diverse physics that one can find 
in real superconductors 7J . Clearly, it is of interest to construct and investigate more general frameworks 
for holographic superconducting systems. Simple generalizations of the HHH model have been proposed 
in [HI ini [To] , by introducing arbitrary couplings and an arbitrary potential which are functions of the 
complex scalar field. A convenient setup is to replace the complex scalar field by a neutral scalar field 
77 and a Stiickelberg field, and then write the most general gauge-invariant Lagrangian which involves 
general functions of the real scalar field 77. In this paper we investigate the main properties of these 
models (reviewed in section 2). We also discuss a further generalization obtained by introducing of a 
generalized "theta term" of the form 8(77) FA F . 

In weakly-coupled superconductors, Anderson's theorem (see e.g. [71 [H]) implies that the thermody- 
namic properties of any superconductor remain unchanged if a perturbation does not break time-reversal 
invariance and does not cause a long-range spatial variation of the order parameter. This basically implies 
that non-magnetic impurities do not affect the thermodynamic properties of a superconductor. There 
are many examples of violation of this theorem in real superconductors, which is itself a sign that the 
underlying physics involves strong coupling. But in general one expects that time-reversal breaking per- 
turbations will give rise to new effects in the system. It is therefore of interest to study mechanisms for 
spontaneous (or explicit) breaking of time-reversal symmetry. Within the context of our general models, 
a natural way to introduce a time-reversal violating perturbation is in terms of a bulk term of the form 
0(77) F AF. Choosing 6(77) such that 0(0) — 0, this term will only be turned on in the condensed phase. 
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While this term does not affect the temperature dependence of the order parameter, we will find that it 
has striking effects in the conductivity. 

This paper is organized as follows. In section 2, we present our general model and derive the equations 
of motion. Section 3 contains analytic results for the behavior of condensate near the phase transition, 
including critical exponents. Section 4 is devoted to the study of the free energy and specific heat near 
the critical point. Section 5 studies the effects of deformations on the conductivity. Section 6 investigates 
Hall conductivities arising from the generalized theta term. A discussion is given in section 7. 

Note added: As this paper was being finalized, [3S] appeared, with some overlap with section 3.2 
concerning analytic expressions for critical exponents. 

2 General holographic superconductors 

We are interested in studying a system undergoing a phase transition where a C/(l) symmetry is spon- 
taneously broken at a certain critical temperature Tc through holographic methods. On general grounds 
[12], the global U{1) current is dual to a gauge field in AdS. Spontaneous symmetry breaking in the 
boundary theory is dual to the Higgs mechanism in the bulk. It is then natural to consider the most 
generic theory under such premise. Extending previous proposals in [HI [9], such generic models have been 
considered in [TO] . They are described by the following U{\) invariant 3 + 1 dimensional^ action 

S = I d'+'xV^g (^R ~ \g{i^) F^'^F,, + _ 1 (a^)2 _ 1 j(^) (9^0 „ . (2.1) 

The Lagrangian describes the dynamics of a metric g^j,, a [/(I) gauge field a real scalar field 77 and a 
Stiickelberg field Q. We work with dimcnsionless fields A^, 77, while \d\ — ~\ (note that [J(77)] = 2). The 
canonical normalization for the kinetic terms can be restored by an appropriate rescaling involving the 
Newton constant Gjsi . The action can also be cast in terms of a complex scalar field 

S = y d^^^x^g (r - \G{i,i,) F^-F^, + ^U{i,i,) - K{i,i,) D'^i^D^^i^ , (2.2) 

where 

^ = f,e^^\ ^ = D, = ^,-^qA,. (2.3) 

The two models p.ip and ()2.2|) are not equivalent in some cases, since ()2.ip permits negative values for 
the scalar field, whereas appearing in the couplings, is positive. In particular, they are equivalent 
when G{i]), U{rf) and J{rj) depend only on even powers of rj with J = q^rj^ + 0{rf ). 

The models typically describe holographic superconducting phase transitions in 2+1 dimensional 
systems. The phase transitions were shown to be either continuous or discontinuous, depending on the 
choice of the G(?7), U{rf) and J{r}) couplings [51 [HI [TO. As we will see later on, the continuous phase 
transitions can be of arbitrary order in Ehrenfest's classification. Interestingly, the string and M-theory 
realizations of holographic superconductors [TH [TJ] belong to this general class of models (12.11) . 

The model (|2.ip admits yet one more generalization (not studied in jlOj) given by 

5S = j d^+^x i 6(77) e'^'^P'^F^^Fp, . (2.4) 

For generic 6(77), this term violates parity and time-reversal symmetry. However, for special couplings of 
the form Q = O^if'^'^^ , with integer fc, parity and time-reversal symmetries still hold if we assume that rj 
is a pseudoscalar transforming as 77 — )■ —77 under P and T . Note that this parity-preserving option cannot 

^This model can be generalized to arbitrary d + 1 dimensions |10) . Here we focus on the case d = 3 for concreteness. 
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be implemented in the formulation (|2.2p in terms of a complex scalar field (in particular, in the model of 
[5]). Similar terms appear from M-theory compactifications, in consistent Kaluza-Klein truncations that 
include a complex scalar field and a real scalar field [13]. In this case, is a specific function of the real 
scalar field. 

For reference purposes, it is useful to define two particular models with couplings having the following 
behavior at small 77: 

Model I : G(77) '^l + mf + 0(77^) , [/(t?) = 1 + ^ r?^ + 0(77^) , J(r;) ^ q^if + 0(77^) . (2.5) 



with 0(77) = 0. The HHH model [5] is the particular case in which k = 0, with U ~ 1 + i 77^ and 
J = q^rj^; the string and M-theory models [131 [HI [H] also belong to this class (see [TO]). It will be shown 
that all these models have mean field critical exponents. Nonetheless, the physics described by models I 
with different couplings will be shown to exhibit strikingly different features. A slight generalization of 
Model I is to allow for odd powers of 77. We will show that including a power 77^^ in the expansion of G, 
U or J already modifies the critical exponents. 

The second class of models that we will consider allows for non-analytic terms in the small 77 expansion: 

Model II: G{v) = I + KT]^ + golvW U (v) = I + + uoM' , Jiv) Q^V^ + jo\v\' ■ (2.6) 

6 

Here a, b, c are assumed to be real, positive numbers, with a, b, c > 2 {G{ri) and J(77) must be positive 
definite for unitarity and uq < for stability). For generic values of a, 6, c, this model contains non- 
analytic interactions in 77. Since the classical solutions we will study have 77 > oH the non-analytic 
dependence on 77 is irrelevant at the classical level. Higher order terms can be added to the expansion. 
We will see that such terms do not modify the dynamics close the phase transition, but can be important 
in other regimes, such as low temperatures. 

A problem of general interest is to determine whether real materials can be approximated by general 
holographic models to some extent. If this is the case, the results in this paper can be used to determine 
the best holographic fit, within our class of models, for a given material. 

Before moving on, we would like to emphasize some important points regarding the interpretation of 
models with non-analytic interactions, such as type II models ()2.6|) in the case of generic (non-integer) 
exponents a, 6, c. It will be shown below that the presence of non-analytic terms in model II allows critical 
exponents to be tuned to arbitrary values in a certain range, thus including values away from the stan- 
dard Landau theory. However, we recall that one does not expect non-mean field theory behavior at large 
Nc- This is consistent with the fact that one does not expect any non-analytic terms in the potential in 
classical Lagrangians originating from string/M theory compactifications, in general. Non-analytic terms 
might effectively be induced by quantum corrections but, since in string theory quantum corrections are 
controlled by 1/iVc, the coefficients of such non-analytic terms would be small in the large Nc approxi- 
mation, vanishing in the classical supergravity limit. Finally, we stress that (as in most approaches to 
holographic superconductors) here we will not be concerned about the UV completeness of the model. 
Type II models (j2.6p with generic exponents may just be regarded as a phenomenological approach that 
can capture some specific features of real superconductors, in particular, to incorporate/parametrize 
certain interesting behaviors close to the phase transition. 

2.1 Ansatz and equations of motion: General setup 

Upon fixing the gauge 6 — the action ()2.ip takes the form 

S = J d'+'x^g R - ^Giv) F^-^F,^ + ^U{ii) ^{dr^f - ^ J(r?)A^A^) . (2.7) 

^In the formulation in 112.21 1 in terms of a complex scalar field the couplings depend on 'ipiji, which is always positive. 
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We now consider the following ansatz 

ds^^-g{r)e-^^''^df + — +r^(dx^ + dy^), A ^ (f>{r)dt , V ^ v{r) ■ (2.8) 
The effective Lagrangian takes the form 

= -2e-^{rgy + ^^G{v)ei <l>'' + ^e-*[/(ry) - ^6"*^'' + ^e* J(7y)0^ (2.9) 
It follows that the equations of motion reduce to 

X' + \rl^ + 2^6^^(77)0^ - , (2.10) 
7 + ^ e"<^'' + ^ + i - i + J(r;)</.2 = , (2.11) 
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= 0, (2.12) 

(2.13) 



In the limit g — > cx) keeping gr/ and qcj) fixed, the matter source drops out from the Einstein equations 
(j2.10p . (|2.1ip and the solution is just the AdS-Schwarzschild black hole, 

3 

5(r)=r2(l-^), x(0=0, (2.14) 

where represents the position of the horizon. In this limit the scalar and Maxwell equations keep the 
same structure, 

._o, (2.15) 



S G J gG{r,) 

^" + ^' + 7) + 2^ ^"^ + T^g ^"^ + ^ = • (2-16) 

In [TU] (generalizing the discussion of [B]) it was found that, even for small values of q, backreaction does 
not alter the thermodynamic properties of the system significantly. Therefore, in what follows we will use 
this no backreaction approximation, with the purpose of simplifying the numeric analysis and obtaining 
some analytical results. 

It is convenient to introduce a new coordinate z = Vh/r, so that the horizon is at z = 1 and the 
boundary is located at z = 0. The equations of motion become 

d^^ + j - _____ ^ = 0, (2.17) 

Oh - O.r, + id.,r + ^^^^ + = ,(2.18) 

where we have defined Lp = (p/rh- In these new variables the asymptotic behavior of rj and is 

, ^^^^ 2 MP -,o^ 

Note that the asymptotic behavior of rj -in turn related to the conformal dimension of the condensing 
operator is due to the special choice of mass in (|2.5p . (|2.6p . which is also the choice of the HHH 
model. The behavior for different masses has been studied in [161 HZl • 
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If 77'^-'^) — ■q^'^^ — 0, then the solution is 77 = in the whole space. As explained in [5J there 
are two schemes that describe spontaneous symmetry breaking. One corresponds to setting rj^^^ = 0, 
such that a global C/(l) symmetry in the boundary field theory is broken spontaneously by a condensate 
(O2) — V^'^^ 0- The other scheme corresponds to 77'^' = 0; in this case the superconducting regime is 
described by a condensate (Oi) = ry^^^ 7^ 0. 

2.2 General models, CFT deformations and quantum critical points 

Let us conclude this section with a few general remarks. We are interested in continuous phase transitions 
where a,U{l) symmetry is spontaneously broken. At the phase transition, correlation lengths diverge and 
the physics is expected to be dominated by a (strongly interacting) conformal fixed point. Universality 
then suggests that few parameters are necessary to distinguish the fixed point. In our case, since ry is 
small near the critical point, we expect the identity of such fixed point to be controlled by the lowest order 
terms in the small condensate expansion of the functions G(ri), U{rj) and Jij]) in (|2.6p . More concretely, 
the dimension and charge of the CFT operators are given by the coefficients of the quadratic terms. In 
addition, as we will see below, critical exponents are dictated by the exponents of the subleading terms. It 
is then clear that our general models have plenty of room for varying the functions without changing this 
data, by changing the coefficients of subleading terms (as we will do in section[5|) or including higher order 
terms. Some of the perturbations, like modifying the potential function U (77) are rather straightforward 
to implement. These higher order corrections can be interpreted as perturbations of the conformal fixed 
point by irrelevant operators. As such, they would not alter the properties of the conformal fixed point, 
but would be in turn relevant away from it when looking at other dynamical properties, among which 
perhaps the most natural one is the conductivity. From this point of view, it is natural to expect that 
these general models might provide a framework for studying quantum critical points, namely second 
order phase transitions at zero temperature as a function of variable couplings|f| 

In general, subleading terms in the G(r]), U(j]) and J{ri) functions represent interactions that might 
effectively account for the result of integrating out massive modes, giving rise to a low-energy effective 
action for rj. From the dual field theory viewpoint, this would translate into an effective free energy 
functional for the low energy effective degrees of freedom, namely the operator (Oi) dual to the rj field. 
This free energy functional will be discussed in section 4. As discussed earlier, non-mean field critical 
exponents arise if one includes non-analytic subleading terms in G{r]), U{rf) and J(?7), and such terms 
should be suppressed in the large Nc limit. Alternatively, one may take a phenomenological approach 
and use the parameters of the model to fit properties of real systems, possibly allowing them to be 0(1). 

3 Behavior near the phase transition: temperature dependence 
of the order parameter and critical temperature 

As shown in [8| IQ} fTO]. the models above exhibit phase transitions at a certain critical temperature Tc- 
Depending on the choice of parameters, these transitions have been shown to be qualitatively different; 
ranging from first order to second order with different critical exponents|f| In the remainder, we will 
concentrate on second order transitions. Previous studies performed in the past have been mainly numeric. 
On the other hand, it would be desirable to have an analytic handle on crucial properties characterizing 
the transition. In this section we describe a simple method to determine the critical temperature and 
an analytic expression of the temperature dependence of the condensate close to the phase transition, 
including critical exponents. 

^ A recent discussion on quantum phase transitions in the holographic context can be found in 1181 . 

*Here in what follows, we use the term second order to refer to continuous phase transitions. In section |4l we discuss 
the behavior of the free energy across the transition in more detail. 



6 



3.1 Critical temperature 

A convenient way to determine the critical temperature of second (or higher)-order phase transitions is 
as fohows. Near the critical point, 77 — > 0, so we can write a perturbative expansion 

V = efi + 0{e^), (^ = ^o + eV + 0(63) , (3.1) 

where fo and fj satisfy the equations 

5,Vo = 0, (3.2) 

To leading order in e, only the quadratic terms in G{ri), U{ri) and J{ri) contribute to the fj equation. 
Thus in what follows we keep only the quadratic terms of model I or model II, and the equation for fj 
becomes linear. The solution for ipQ is simply 

^o = ^{l-z) = ^i---). (3.4) 



Substituting the expansion p. 51) into (13.31) . we find that the differential equation p.3p depends only the 
combinations k = n/q^ and A = ^q/r^. The boundary conditions for the 77 equations are 



fl{l) = rih. f,'{l) = \r,u^ {K~q^)^\i^h^' {k-l) . (3.5) 



Since the rj equation is linear, the solution is of the form 

f]^7^hF{A,k; z) . (3.6) 

Expanding near z = 0, we have 

^ ^ci{A,k) z + c2{A,k) z"^ + .... (3.7) 

In the scheme where (O2) is non- vanishing, the critical temperature is then determined from the require- 
ment that 

ci(^,k)=0. (3.8) 
In the scheme where (Oi) is non- vanishing, the requirement is 

C2(^,k)=0. (3.9) 

In either scheme, these conditions determine A = A{k). Then one uses T — 3r^/47r and r/j = ^-q/A or 
rh = yfpqJA. This gives a critical temperature of the general form 

Tc^^ qt^f{n/q^) , (3.10) 
47r 

or 

3 



Tc = ^V^P ViW?) , (3.11) 



which are appropriate for calculations at fixed ji or fixed p, respectively. Here / = 1/ A{k). 

Since the method implies solving an ordinary second order differential equation, the critical temper- 
ature can be determined with high accuracy. Here we use the scheme where (O2) is non- vanishing. In 
particular, for fixed ^ = 1, wc find the value Tc{k. = 0) = 0.05874734 q (or Tc[k = 0) = 0.1184267 ^ if 
we work at fix charge density p = 1, in agreement with [5]). A plot of as a function of n/ q^ at fixed q 
is shown in figure [TJ We see that at large k, T'^ becomes linear with k and, for large negative k, tends 
to zero monotonically. 

Thus, the parameter k can be used to tune the critical temperature to any desired value. A similar role 
is played by the parameter q or the scalar (which for our models was set to the special value = —2), 
i.e. they also affect the critical temperature. 
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Fi gure 1: T'^ vs. k in second (or higher) order phase transitions in Model I or Model II ((/ — /i — 1). 
3.2 Critical curve for the order parameter 

As numerically shown in [HI |9l |10| , the model exhibits second order phase transitions in certain regions of 
its parameter space. In second order phase transitions the order parameter continuously approaches zero 
as we approach the critical temperature from below. Translated into our holographic approach, close to 
the phase transition the scalar rj should be arbitrarily small. A useful approach for studying the vicinity 
of the phase transition is then to solve the equations of motion (I2.17P and (|2.18p in terms of a series 
expansion of j] and Lp near the horizon, located at 2 = lH 

ryW(z) = 77^ + 771(1 -z) + 7^2(1 - z)2 + ... + ??Ar(l - z)"^ , 

(^W(z) = ^i(l-z) + ^2(l-^)' + ... + ^iv(l-z)^ , (3.12) 

where we truncate the series at some order N . 

The coefficients 77^ and (pj can be solved in terms of rjh and (pi using the equations of motion for every 
j > 1. The series p.l2p converges to the solutions ri{z) and (p{z) as iV — > 00. Then, we match these 
expansions with the asymptotic behavior at the boundary (|2.19p 

77^(0)^7700(0) r ^w(o) = (^oo(o) 

{ (3.13) 

9,77^(0)^9,7700(0) [ 9,^W(o) = a,^oo(o) 

We consider model II, which has eight parameters: k, q, go, uq, jo, a, b and c. As we will see, model 
II captures all the terms that are important close to the phase transition for a completely generic model. 
For concreteness, we focus on the (O2) = scheme. There are five parameters rih,(pi (from the horizon 
side) and 77^^^ /i, p (bom the boundary side). These conditions can be solved for four parameters in 
terms of either p or Furthermore, close to the critical temperature, Oi — t?*-^-* is linear in rjh- 

Expressing in terms of T using — AirT/S, we find the following generic structure 

1 - = AN{K,q^){Oif + goBN{a){Oir-^ + uoCN{b){Oi)''-^ + joDN{c){Oiy-^ + . . .(3.14) 

where ". . ." denotes terms that are of higher order in 77/1. The key point is that the structure of p.l4p 
remains the same at each successive order in N. The different functions (including Ti^-*) get corrections, 

^An alternative method that leads to the same result follows the idea of the previous subsection, namely to look in the 
vicinity of Tc and solve the differential equations in perturbation theory in powers of rj to the next order. 

^Another possibility is to implement the matching at some intermediate point Zm, with < Zm < 1, as done in 1191 . 
While having variable Zm allows one to obtain better solutions at fixed order A' in the expansion, it introduces complicated 
algebraic equations that prevent the application of this method at large Nc- 
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but the general form of the equation, which remains unchanged, aheady reveals interesting information on 
the phase transition. In particular, it gives the explicit functional dependence of the order parameter on 
the temperature, hence it gives the explicit analytic expression for the critical exponent (3. For a, b, c> 4, 
the leading term in p.l4p is jf. In this case we find 

(Oi) = const, (l - f ) ' ' for T = , (3.15) 

as in mean field theory. When either a, 5 or c is less than 4, then the leading term is 77°""^, where 
ao = Min{a, 5, c}. In this case, we find 

f tV 1 

(Oi) ^ const. 1 - T^r , P = ^ , ao - Min{a, b, c] . (3.16) 

Because of the presence of the term AM{K,,q^){Oi)'^ , generically one has the bound /3 > 1/2. This 
value, as well as the possibility of deriving it from an expansion like p.l2p . was already advanced in 
[SJ 110) and is in perfect agreement with the numeric results. Remarkably, the correct value of exponents 
can be determined from the expansion at lowest order, in contrast to quantities such as the critical 
temperature. Their values are robust provided the corresponding coefficient in (j3.14l) does not go to zero 
in the N ^ 00 limit, something that can be easily checked. Recently, the same value of the critical 
exponent (3 was independently derived in [55], where it was discussed in the context of general U{rf), 
J{r]). Furthermore, in the specific case in which the Breitenlohner-Freedman (BF) bound is saturated, 
|28| found an analytic solution from which the numerical coefficients in p.l4p can be determined. In 
sectional we explain how the same result can be deduced from minimizing the free energy. 

In principle, it is possible that An^oo{k, q^) vanishes for some real value of n/q^. If this is the case, 
one might think that in this special point of parameter space one can obtain models with /3 < 1/2 by 
taking a,b,c > 4. However, numerical study suggests that before this occurs the transition changes into 
a first order phase transition. Finally, let us mention that an equation with the same structure as p.l4p 
is obtained in the rj^-^"^ = scheme where (O2) takes an expectation value, so the same considerations 
apply in this case. 

An important feature we learn from the critical curve p.l4|) is that Tc depends only on n and g, 
i.e. on the quadratic terms in the expansion of G{ri), U{'if) and J(?7). By introducing a general mass 
m? parameter in U , one sees that Tc is more generally a function of k, to, and q (the mass dependence 
of Tc was previously discussed in detail in |17)). This important feature is confirmed by the analysis of 
the previous subsection, where we saw that the critical temperature is determined by a linear differential 
equation and it is thus unaffected by higher order terms in Girf), U{ri) and J{'q). In particular, this 
feature was the content of conjecture 2 of [TO] . 



4 Free energy 

Let us now turn to the study of the free energy for our general models. As standard, the free energy is 
given by the gravity on-shell action. This quantity is divergent, and thus needs to be regularized before 
further physical interpretation. Regularizing the action with a radial cut-off r^, it is not difficult to show 
that the structure of such action is 

+ + (4.1) 

where the dots stand for finite terms irrelevant at this stage. In addition, we keep the superscripts to 
remind that we will focus on the case with fixed chemical potential and the (O2) = scheme. Since this 
action is divergent, we need counterterms to holographically renormalize it. Given the structure of the 
on-shell action, it is nevertheless clear that the same counterterm as in [8] will do the job 
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ASf 



(4.2) 



being 7 the induced metric on the boundary B. Then, the free energy for the general models (|2.ip takes 
the form 



2 2 



^ / dz 



V 



l-{r^d,U-2U), 



2""'" ■ 2z2(l~z3) z4 

Near the critical point, where rj is small, we can use the expansion for type II models 



(4.3) 



2 



Oi O2 



'-fid. 



(1 



and write 



where we have introduced x — ^ ^ V dropped the superscript. We have also set k = to maintain 
the same critical temperature. 

In order to further proceed, we can use the series solution to the equations of motion as in p.l2p . 
Inserting the resulting expressions in (|4.4p we obtain 



M/W = W^''\T){n,q^T) + w!,''\n,q\T)Ol+g,W(^\n,q\T,a)Ol 
+ uo Wi""^ {k, q', T, h) 0\ + JO W^i"^' {n, q\T,c)0'i + ... 



(4.5) 



where ". . ." denotes terms with higher powers of 77^- We keep this schematic form since, once again, the 
particular expressions of the coefficients VK/^'' are cumbersome and not too illuminating. As for p.l4p . 
the general structure of (|4.4p is independent of N . Indeed, an identical discussion applies also here. We 
have also verified that the expansion in (|4.5p approximates the exact (numerical) free energy close to the 
phase transition. 

The critical temperature at order N corresponds to the value at which W^2^^ changes its sign 

V7i^)^(T-rW(«:,g^)) + ... (4.6) 



(«,g2). 



Wq^-* and Wc^'' remain finite at Ti^-' 

As discussed in section 2.2, the holographic free energy functional (|4.4I) can be thought of as a sort of 
generalized version of the Landau-Ginzburg free energy. For integer values of the exponents a, h and c, we 
can think that the extra couplings in G{ri), U{rj) and J{rj) correspond to the effect of higher dimensional 
operators in the bulk effective field theory. Such terms are known to arise in consistent truncations of 
String/M-theory [141 115) . Further exploring their possible origin in the boundary field theory is of great 
interest, but beyond the scope of this paper. 



4.1 Specific iieat and higher-order phase transitions 

The free energy (|4.4p contains all the relevant information about the system close to the phase transition. 
By minimizing it with respect of Oi we can obtain the physical VEV of the condensate as a function of 
the temperature. In particular, we can consistently check that the critical exponent /3 obtained in this 
way coincides with the one we got in the previous section. Substituting such value into the free energy 
(|4.5p . we can determine the behavior of the specific heat close to the critical temperature. We obtain 

= -^?^^ const. (T, - T)-" , a--i^. (4.7) 

where oq = Min{a, h, c} appeared earlier in (I3.16p . Note in particular that we automatically get 
a+2/3 = 1. This result was also recently derived in |28| for the case with non-trivial U{rf} and J(??), simply 
using the scaling relations. Our result holds more generally and, since our derivation is independent of 



10 



-17 -16 -15 -14 -13 

T 

Log[l ] 

Figure 2: Ac„ as a function of temperature close to the phase transition for J = 7]^ — i]'^ + rf, with 
c = 3 + /c/5, fc = 0,...,4. 



the scahng relations, it provides an independent proof that they are satisfied as we discuss below. Figure 
[2] shows Ac„ as a function of temperature close to the phase transition for non-trivial J{rj). 

In I^D], it was shown that holographic superconductors have a universal value 7 = 1 for the critical 
exponent associated with the thermodynamic susceptibility, xt (T^ — T)~'^ . The same analysis leads 
to 7 = 1 also in the present more general models. Putting these values together, we find that the critical 
exponents of the general Model II verify the Rushbrooke identity: 

a + 2/3 + 7 = 2. (4.8) 

Since 2 < oq < 4, we find that /3 > 1/2 and a < 0. For a < (or 2 < oq < 4) the specific heat is 
continuous across the transition. This means that the transition is at least of third order (using the 
Ehrenfest convention that the order of the phase transition is the lowest one at which the derivative of 
the free energy is discontinuous). More generally, the phase transition is of order n for 

n - 1 > 2;3 > n - 2 . (4.9) 

In other words, for these models the order of the phase transition is given by 

n=r2/3 + ll = r^l , (4.10) 

where \x\ denotes the smallest integer greater or equal than x. The transition will be of higher order 
than 3 if ao < 3. However, in this case, /3 > 1 and the critical curve becomes concave near T^. These 
cases are of some theoretical interest because there is a temperature Ti < where d?{Oi) /dT"^ changes 
sign. The change in the sign of d?{Oi) /dT^ could be an indicator of another phase transition, though we 
will not further inquire on this point in this work. 

It should be mentioned that there are examples of real superconductors with third order phase transi- 
tions (see e.g. [21[). We are not aware of any real superconducting system with well established transitions 
greater than third order (the experimental data given in [21] are also consistent with the assumption that 
the transition is fourth (or even higher) order; distinguishing the two cases would require high exper- 
imental precision data which to our knowledge is not currently available). In our context, third order 
phase transitions can be described by simple models of type II (j2.6p with the addition of rf interactions 
in G{vi), U{ri) or J{ri). 

For ao = 4, we have /? = 1/2, a = and there is indeed a jump in the specific heat, which can be 
compared with the BCS result. For concreteness, let us focus on the case in which only J — rj^ + jq 
is not trivial. Computing the jump in the specific heat, we obtain 

AN) „-2„4 , „(^V) ■ „8 , „W „12 

Ac.(r.) = c.(T.) - c„(T.) = i]'-'''^ ' , (4.11) 

(r -a^ jor 
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where the a,- are numerical factors whose values depend on the order N to which we approximate the 
solution. The functional dependence of the specific heat on q and jo does not depend on N. In particular, 
at jo = 0, i.e. for the HHH case, it is a constant independent of q. Furthermore, consistently with the 
analysis in section 3.1, we also find that Tc ^ q, in fact Tc/q ~ 0.0587 (see (I3.10|) and below). Using 
numerical results to fit the proportionality coefficient, we obtain 

Ac,«— T,. (4.12) 
Interestingly, in BCS theory, the analogous quantity evaluates to 

Ac, = ^^7V(o)Te«9.38 7V(0)r,, (4.13) 

where A^(0) is the density of states at the Fermi surface. Comparing with (4.12) we see that, in this 
particular model, 1/q seems to play an analogous role as A^(0) (although there is of course no Fermi 
surface here). 



5 Conductivity 



So far we have been mainly concerned with the universality properties of the phase transitions described 
by our generalized models. As noted in subsection (2.2), these are controlled by the lowest order terms 
in the small condensate expansion of the functions G{r]), U{r]) and J{r]). Higher order terms can be then 
thought of as irrelevant perturbations of the conformal fixed point, thus not modifying the universality 
properties of the phase transition. However, they do play a key role in dynamical aspects away from 
Tc- Indeed, our models have a priori enough room to accommodate a plethora of behaviors that might 
model real systems. One of the most natural such dynamical quantities to consider is the conductivity. 
By the Kubo formula, it is related to the spectral density, which in turn represents the density of energy 
eigenstates at energy w, weighted by the overlap with the electric current operators. 

Let us first briefly review how to study the conductivity in our general models. Following [4l[5l|6], we 
consider time-dependent perturbations of = ax{r) e~^^* and gtx = f{,r) e~*"*. These fluctuations are 
governed by the following equations 



+ — 



2 



G 



_J_ 

'gG 



f'-- f + G cj^'a, = 



9 



f' + -f 
r 



(5.1) 
(5.2) 



Substituting the second into the first equation, we find 




GcT 
9 



9G, 



(5.3) 



The asymptotic behavior of the perturbations is found to be 

/(I) 



(5.4) 



The conductivity can then be obtained from the formula 



Jx 



.(1) 



,(0) ' 



(5.5) 



where in the second equality we have used the AdS/CFT dictionary. 
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Following the ideas in section [^?^ we now investigate in detail the behavior of the conductivity under 
deformations that preserve dimension and charge of the condensate and the critical exponents. A simple 
implementation of this idea is given by the choice: 



G = 1 



6 



J = rf + joT] , 



(5.6) 



with Jo > 0. This choice provides a one-parameter family of deformations where J remains positive for 
all rj. For any given temperature, we can get arbitrarily close to the HHH model by reducing Jq. Figure 
[3] shows the conductivity for jo = 0.6 in the (Oi) = scheme (in this section we work at fixed charge 
density p) . The most salient feature of this plot is the appearance of resonance peaks in the conductivity, 
which increase in number and become narrower and higher as the temperature is gradually lowered. A 
similar effect is also seen at fixed temperature by increasing jo. One can also see that increasing jo leads 
to an increasing of the height of the peaks. The lowest temperature curve in fig. 3(a) is at T = 0.24. Here 
the first peak in Re((T) cannot be seen numerically because it has become narrower than the numerical 
grid. However, its presence can be inferred from the l/(w — uj') behavior of Im(cr), shown in fig. 3(b), 
through the Kramers-Kronig relation. 

Consider now two different deformations: 



gov 



6 



J = 77% 



(5.7) 



G = 



1 



U 



1 + 1 



(5.8) 



1 + 9oV^ ' 

with go > 0. The three models above belong to the same universality class, since they have the same 
critical temperature and the same critical exponents as in the HHH model (obtained by setting jo = go = 
0). A natural question is to what extent transport properties are sensitive to the choice of deformations. 
Figure S] compares the conductivities of the three different deformations (15. 6p , (|5.7I) and (15. 8p in the 
(Oi) = scheme at a given temperature T = 0.415 Tc- One can see that the conductivity undergoes 
significant changes relative to the HHH model, with strong dependence on the specific deformation. At 
this temperature, the model (15. Sp exhibits sharp peaks while the models (15.61) . (j5.7p have a smoother 
behavior, closer to BCS. In particular, the figure shows that a small deformation of the model can change 
the density of energy eigenstates of the system in a dramatic way. 
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Figure 3: (a) Real part of the conductivity as a function of frequency for J = rj"^ +jov'^ with jo = 0.6 in the 
(Oi) = scheme. The curves correspond to different values of T/Tc equal to 0.24, 0.29, 0.50, 0.61, 0.81 
(the curves with lower temperatures are those that go to zero more rapidly as w — > 0). (b) Real and 
imaginary part of the conductivity for the same model at T — 0.20. 
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Figure 4: Conductivity as a function of frequency for HHH (dashed-dotted line) and for the different 
deformations given by the models (|5.6p . (15.71) and (15.81) in the (Oi) = scheme at T = 0.415 Tc- 



Similar resonance peaks have been observed in [TB] by varying the dimension of the condensate (i.e. 
the mass of the scalar field), making it approach the BF bound. Our results show that it is not necessary 
to set the mass of the scalar field to any particular value. The same phenomenon can be obtained 
by higher order modifications of the functions G{rf), U{r]) and J{ri). In fact, this is not surprising; the 
conductivity peaks are a low temperature effect. As such, they are controlled by the value of the functions 
for large values of the condensate in contrast to e.g. critical exponents, which are determined by the small 
condensate expansion. In the simple models of |16| . which fit into our class I of models, tuning the mass 
of the scalar m (a rather drastic modification of the dual CFT) is the only available way of controlling 
the behavior of J for any value of the condensate. 

From the bulk perspective, the mechanism that gives origin to this peaks was clarified in [23]. We 
now generalize this argument to our context. Introducing a new coordinate 

eX/2 

du = rfr, (5.9) 

9 

and defining = \/G ax , ()5.3p takes the form of a Schrodinger equation 

-^ + V{u)^ = uj^-9 , (5.10) 



with 



where prime denotes derivative with respect to r. It should be noted that V not only depends on the 
couplings G, J, but also implicitly depends on U, through the dependence on rj, (which are in turn 
determined by the coupled system (|2.15|) . (|2.16|) '). Note also that in the present approximation, where 
back reaction is neglected, x can be set to zero and the first term proportional to can be dropped 
from (fSHj) . 

Now consider for example the model (|5.6p . The potential V is shown in figure [5{a), for jo ~ 0.05, 
along with the HHII case, jo — 0, always in the scheme where (O2) 7^ 0. 

The conductivity is obtained by sending an incoming wave from the right, which will be partly 
transmitted and partly reflected, with a reflection coefficient TZ. As observed in [23] the conductivity 
(j5.5p is nothing but 

1-7?. 

a(w) = . (5.12) 
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Figure 5: (a) Potential for the model (|5.6p with jo = (dashed line), corresponding to the HHH model, 
and for jo = 0.05 (solid line) in the scheme where (O2) 7^ 0. The temperatures are 0.78, 0.24, from 
bottom to top. (b) Maximum of the potential as a function of the temperature. The dashed line is the 
model (I5.8P with G = 1/(1 + Q.lrf)] the solid lines correspond to (from bottom to top) HHH model, 
J = ff + QMrj'' and J ^ if + 0.2r]^. 

If the potential is very high compared to (and as long as the area under the potential is also large), 
the absolute value of the reflection coefficient will be close to 1. The peaks are then produced when the 
reflected wave has a relative phase equal to {2n + 1)tt, so that 1 + 7?. is near zero. Figure [5] (b) shows 
the maximum of the potential Knax as a function of the temperature for different models. We see that 
Vmax increases as the temperature is decreased. This gives room to new resonance frequencies satisfying 
the condition lu <C y/Vmax and hence the emergence of new peaks. In HHH, in the scheme where (O2) is 
non-zero, there is a small region where the incoming wave is oscillating near u — and as a result there is 
a small enhancement in Re(cr) at a certain frequency. Furthermore, Vmax tends to a finite value at T = 
which is not big enough to allow for a reflected wave with \TZ\ ^ 1 and thus for the formation of peaks. 

It is also worth noting that, as long as the potential stays finite at all temperatures (including T = 0), 
there will be no hard gap in these models, since the probability of tunneling will never be zero. As a result 
\TZ\ < 1 and the conductivity will have a finite (albeit exponentially small) value at low temperatures. 
However, figure [5] (b) suggests that for some models Knax might actually go to infinity as T — ^ 0. If 
this is the case, there could be a hard gap and an infinite number of resonance peaks in this limit. The 
figure [5] (b) also explains why the model G = 1/(1 + 0.1?7'') has a more suppressed conductivity at low 
frequencies and narrower peaks at a given temperature, see fig. HI as we can see, in this case, at a given 
temperature the potential barrier is higher. We will make further comments on the issue of peaks in the 
Discussion section. 

6 Models with dj^Q ^ and Hall effect without magnetic field 

We shall now incorporate the term (|2.4p . Even though a priori it contributes to the equations of motion 
as soon as 9,,0 ^ 0, it is easy to see that for the particular ansatz (|2.8p it actually gives a vanishing 
contribution. Therefore we find exactly the same uncondensed and condensed black hole solutions irre- 
spective of the coupling 8(77)0 Nonetheless, this coupling will affect the conductivities in an important 
way. The equation of motion for Ay shows that this new interaction leads to the interesting effect that an 
electric field in the x direction turns on an electric field in the direction y, in very much the same fashion 
as in the Hall effect, this time without the presence of a magnetic field. The basic idea behind this effect 
was first raised in [2^ in the context of a model having a constant Q. In that case the Hall conductivity 
is a numerical constant, present even at zero temperature and zero charge density. The advantage of a 

^The dyonic solution is also unaffected by the © term in the uncondensed phase, but it changes in the condensed phase 
if dr,0 + 0. 
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coupling &{ri) with 9^0 7^ is that it can incorporate non-trivial temperature and frequency dependence 
on the Hall conductivity a^y through the non-trivial coupling to the scalar field. Moreover, we have the 
option of choosing 8 such that 0(0) = 0; in this case the Q term will be turned on only in the condensed 
phase. More recently, in [T3], a non-trivial Hall effect was found by essentially the same mechanism as 
in this paper, arising from a term of the form f{h) F A F, where h is an additional real scalar field and 
/ is a specific function emerging from the compactification. A different way to obtain Hall conductivity 
appeared in the context of p-wave superconductors [25] . 

6.1 Hall effect in the general model 

Let us consider the most general model with non-trivial G{ri), U{ri) and Jijj). The relevant part of the 
action containing the U{1) gauge field: 



S = Jd^x (^v^ {-\Giv) F^^'^F^, - i J(77)A^A^^ + ie(77) e>^'^P^F^,Fp}j . (6.1) 

The variation with respect to leads to a generalized version of the London equation 

dt^^g G{r^)F'^- - J(r,) A'^ - e^'^P'^ F^,dpQ{fi) = . (6.2) 



Using the ansatz (|2.8I) we find that the equation (I2.12p for the potential (p remains unchanged. In order 
to calculate the conductivity we consider a perturbation = A^(i, r) to the background. We will ignore 
backreaction and consider as usual the linear order. The effect produced by the new O term in ()6.2|) is to 
couple (at first order) the bulk perturbations A^ and Ay. Indeed the equations of motion are (e*''^*' — 1) 



d^A^ + 1 ^ + ^ \drA,- -dfA^ - ^ A, - ——dtAydrQ = , (6.3) 
\ g G J gG gG[ri) 

dlAy + ( ^ ^ - ^^dlAy -^Ay + -^9,^,9,9 = . (6.4) 



Expanding in Fourier modes 



A^= I duj e~"^^a^{r;uj) , Ay = dio e-*'^*ay(r; w) , (6.5) 



the equations for the modes ax{r;uj)^ ay{r\uj) become 



G' \ , /CJ2 J 0/ 



~2 7^< I ^ 2^0-1 = , (6.6) 



g • G r-' \g^ gGn^^'^gGir,) 

g' G' X , fuj^ J \ 
- + ^ a. + ^ Tz]a„- 



< + ( - + - j < + - - ja, - 2..-^ a. . , (6.7) 

where prime indicates differentiation with respect to r (so that G' = djjG rj' , 8' = 9^0 rj'). The system 
p.6p - (|6.7p can be decoupled introducing complex coordinates, z ^ x + iy, so that 

Az = ^(a^ - iay) , Ag = ^{a^ + iay) . 

Using the notation, 

the equations take the form 

PA, - QA, = , PA, + QAj ^ . (6.9) 
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Note that the asymptotic behavior of (|6.9I) is not modified by the new term; indeed for r — >■ oo we have 
A'^ + - 26.^77' A, = , + -^A', + 2uj^ii A,^Q. (6.10) 

Since 77' is 0(l/r^) or 0(l/r'^), the 8 term can be neglected at large r. Then the asymptotic solutions 
are the same as in the = case, 

A, = ^(0) + , A, = Af + ^ . (6.11) 
r r 

Finally we observe the fact that the equations (|6.9p do not depend on the value of 8(0). 

To compute causal behavior, we solve for the fluctuations with ingoing-wave boundary conditions at 
the horizon. This requires 



1-— ) a,{r) , A,^ C-Jl~—) a,{r) , (6.12) 



with az{r) = 1 + a^z\^ — r/vh) + . . . and similarly for Because (j6.9p are linear and homogeneous 
equations we can first set Cz — Cg — ^ and find two solutions Az^ Ag- Then the most general solutions 
are obtained as 

Az = CzAz, A-z = C-zA-z . (6.13) 

From the above expressions we find, 

= C-zA-z + CzAz , a„ = {C-zA-z - CzAz) ■ (6.14) 
To uncover the physical meaning of Cz and Cg, we consider the asymptotic behavior (|6.1ip and find 

«.(..) = 4") + ^ , a,(a.)=aW + ^, (6.15) 

4") - C-zAT + CzA^^ , a(i) = C,Af + , (6.16) 

a(°) = -^{C-zAf - C,4°0 > <^ = -i{C-zAP - . (6.17) 

As is standard in AdS/CFT, the leading term in the asymptotic expansion of the fields is related to the 
source in the dual theory; therefore, from the definition Ei — —dtAi, i = x,y, we obtain the system, 

= iwai"^ = CziujAf^ + CziujA^°^ 



(6.18) 

) 



iE^ = i^ujai°'> = CziiuAf'' - CziojA^°'' 



In this way, the integration constants Cz and Cg get related to the physical sources E^ and Ey. Solving 
for Cz and Cg we find the expressions 

Cg = , Cz = . (6.19) 

The asymptotic coefficients ai^"* and a^"^ in the expansion (|6.17p can now be written in terms of the 
electric field components as 

' 1 4°' 4°V 2 " \Af> Af j 2 ■ ' ' 



^Equations l|6.9p contain two different sources of singularities at tfie horizon, coming from inverse powers of ~ 
(1 — r/r/i)'^ and g ~ (1 — r/r^,). The exponent —iu)/3 takes care of the first one, while the second one is canceled by 
choosing appropriates values for ai^^ and a^^\ This determines one integration constant in l|6.9| l. whereas Cz and Cz 
remain as free parameters. 
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,(1) 



l(i) 



.(0) .(0) 

^ J-Vy. -^Z 



.(0) .(0) 
■^z 



(6.21) 



This expression will be the starting point for the discussion about the conductivity. 
It is convenient to rewrite the action (16.11) in the form 



S= d 



1 



1 



1 



From (16.21) we have the relation 

^Jiv)A^A^ 



^d^{^G{ri)Ff^'')A, + F^.dpQA, = 



Substituting into (|0^ we find 

1 



So.s. — ' 



where So.s denotes the on-shell action. We now specify our ansatz. We obtain 

1 



(6.22) 



(6.23) 



(6.24) 



So.s. 



- 1 d^. 



g G{^)F'--A, - Qij^yp-F^pAo 



g{r){A^drA^ + AydrAy) + Q{ri){AydtA^ - A^dtAy) 



(6.25) 



Using the boundary condition Ax{rh) = Ay{rh) = and the Fourier representation ()6.5p we find 



So.s 



-— I d X I dcj 



g{r){ax{r,Lo)drax{r, -w) + ay{r,uj)dray{r,-uj)) 



d^x I duj Q{rf) iujay{r,uj)ax{r,—uj) — iujay{r,—uj)ax{r,uj) 



1 
2 

Finally, substituting the asymptotic behavior for the fields, we obtain the result, 
1 



So.s. 



d X dcj 



a(")(a.)a(i)(-c.) +4°) (c.)aW(-..) +^^e(0)e»^af)(c.)af (-^) 



(6.26) 



(6.27) 



where i,j — x,y and we use the convention e^^ — 1. 

In the AdS/CFT dictionary, the leading term in the asymptotic expansion of the fields determines 
a source in the dual theory, while the "normalizable" term will give the expectation value of the dual 
current. Therefore, from (|6.17p and the definition Ei = —dtAi, we get. 



Jy — 



SSo.s. 

oa'x 

SSo.s. 



5a. 



= a(i)(a.) - za;e(0)af (..) = a^^\u:) - Q{0)Ey 
aW(c.) + zc.e(0)4o)(c.) = a«(L.) + Q{0)Ex 



(0) "'y 



(6.28) 
(6.29) 



Using the expressions (|6.20p - (|6.2ip we derive the components of the conductivity matrix, defined by 



Ji — '^ij Fj , 



2i 



Wzz — CT22) 



(6.30) 
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where Uzz, (^sz represent the conductivities for left-oriented and right-oriented circular polarizations of 
the electric field, 

1 A^z^ 1 

(^zz = -. 77^ , a^z = ■ (6.31) 



y-vz "^z 

The relations a^x — <^yy and a^y = —<^yx are a consequence of isotropy. In a parity-preserving theory, a^z 
and Ozz are equal to each other, and the Hall conductivity a-j^y (proportional to the difference) vanishes. 
On the contrary, in a parity-violating theory, they differ, giving rise to a non-trivial Hall effect: turning 
on an external electric field in the x direction implies that the system must automatically produce an 
electric field in the y direction and a non-trivial Jy current (and viceversa). For a generic 6(?7), parity 
symmetry is explicitly broken in the present model and we indeed obtain a non-trivial Hall conductivity 
axy^ even if 0(0) = 0. In turn, as described above, for special couplings of the form O = Q^if^^^ , with 
integer fc, parity and time-reversal symmetries still hold if we assume that 77 is a pseudoscalar transforming 
as 77 — >■ —r\ under P and T. However, in this case, in the condensed phase the non-trivial profile of 77 
spontaneously breaks the parity and time-reversal symmetries, thus leading again to non-trivial Hall 
conductivity. 

In the uncondensed phase, r\ = 0, and the parity-violating interaction is not turned on. The system 
has vanishing Hall conductivity as expected. 

6.2 Numerical analysis of the conductivities 

Here we shall consider as an example the HHH action '5' with the addition of the term (|2.4I) , with = 77" . 
The total action reads 

1 f... l".-v^^+ 6 ^^^1 ^^^^,_1 _^ ^^^^ ^^^^ 



So = T^y^^^^-(^(^-i^^^^.^ + Il(l+6^^)-2(^'?)^-2'^V(5.^-^^ 



77" e^'^P'^F^^Fp, ] . (6.32) 



We take q = 1 and L = 1 and, once more, fix the gauge 9 = 0. The n = 1 case is of interest, being the 
simplest case, and also because P and T symmetries are preserved. The n = 2 case is also of interest, 
since this case is easily incorporated in the context of HHH model with a complex scalar field i/j by adding 
to the HHH model the interaction F AF. 

Our aim is to compute axx and axy using (j6.30p and (I6.3ip . by numerically solving the differential 
equations (|6.9p . which for this particular choice of couplings read 





(- 




- 271W77"- 


-1^ 


9 




9 




9 




(- 






-1^ 






9 




9 



-z=Q . (6.33) 

with g given in (|2.14p . Solving these equations numerically with the boundary conditions as described 
in the previous section, we compute the conductivities (|6.30l) and (|6.3ip at different temperatures. The 
results for the model with 71 = 1 are shown in figures [nKa),(b), obtained in the scheme where (Oi) = 0, and 
in figures [7Ua),(b), obtained in the scheme where (O2) = 0. At high frequencies, ^e{<Jxx) approach the 
normal phase behavior seen at T = Tc, showing that the mechanism for conduction is through the normal 
phase charge carriers, while Ke{axy) approaches zero in this high frequency regime. As the temperature 
is lowered, there is a gap in the direct conductivity Re{(7xx) B,t lu < ujg, with ujg ^ ^ (O2) and ujg ^ (Oi) 
in the two different schemes. This is expected, lOg representing the minimal energy to break the Cooper 
pairs. On the other hand Re(tTa;j,) has no gap; it approaches a finite value as w — >■ 0. In the limit that the 
temperature goes to zero the curve V>js,{axx) injHl^a) exhibits a pronounced resonance peak at a frequency 
a; slightly higher than Wg. We also see a mild minimum and a mild maximum at some higher values of 
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the frequency. In the other scheme where (Oi) 7^ 0, wc find a smooth behavior after the gap. The model 
with n = 2 reproduces qualitatively similar features but with many more peaks. 

We have also verified the following important feature. While the direct conductivity Re(o'a;a;) has 
the expected delta function singularity at a; = 0, the Hall conductivity Tie{axy) has a finite value in the 
DC (w = 0) case. This is seen more clearly from the behavior of lm{axx) and lm{axy) near w = 0: by 
the Kramers-Kronig relation, a pole in the imaginary part of the conductivities implies a delta function 
singularity in the real part (see e.g. [5]). The numerical analysis shows that lm{ijxx) ~' 1/w while 



Im((Ta;y) 



' w as uj 



0. 




(a) 



(b) 



Figure 6: (a) Yle{axx) vs. uj (lower temperatures corresponds to curves with lower intercepts at w = 0) 
and (b) Re(cra;y) vs. w for the model (|6.32p with n = 1 at various temperatures, in the scheme where 
(O2) is non-zero (here the curves with lower temperatures are those which have higher value of 'RBia'xy) 
at the frequency of the peak). 




(a) 



(b) 



Figure 7: (a) ^e{axx) vs. u and (b) ^e{uxy) vs. u for the model (j6.32p with n = 1 at various temperatures, 
in the scheme where (Oi) is non-zero (same conventions as in figure [6]). 

Having a finite DC Hall conductivity, it is of interest to investigate its dependence with the temper- 
ature. This is shown in figure |8] for the n = 1,2,3,4 models. The first obvious feature is that the Hall 
conductivity vanishes at T > Tc, since in this case r] — and the generalized "theta term" does not 
contribute. Near Tc, in general the DC Hall conductivity has the behavior 



Re{axy){uj = Q)^{T,-Ty 



(6.34) 
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where v depends on n (e.g. v ~ 0.4 and v ^ \ for 71 = 1,2 respectively, and v > \ for n > 2). In 
the n = 1,2,3 cases, at lower temperatures the DC Hall conductivity is approximately temperature 
independent. 
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Figure 8: ^e.{ij-j;y){yi = 0) vs. T for the model (|6.32p with n = l, n = 2, n = 3 and n = 4 (from top to 
bottom). Since the different curves correspond to different theories, for comparison Re((T2;i,)(w — 0) has 
been normalized to 1 at the lowest temperature numerically available. 



7 Discussion 

In the first part of this paper we have explored the properties of the holographic superconductor system 
described by the action ()2.1|) . for different couplings G(?7), U{rj) and J{rj) with the small condensate 
behavior ()2.5|) or p.6|) . We confirmed that models of the type I (|2.5p have second order phase transition 
with mean field theory exponents, while models of type II can have phase transitions of arbitrary order 
and critical exponents that depend on the parameters of the model. We analytically determined this 
dependence finding the following behavior 



{0)^{T,-Tf, Ac, -(Te-T)-" , ^^^(T.-T)-^, (7.1) 

1 4-ao 

, a = 

ao — 2 ao — 2 



with 

an = Min{a, 6, c} , 2 < < 4 . (7.3) 

It follows that the standard Rushbrooke identity among critical exponents, a + 2/3 + 7 = 2, is satisfied. 
This is a reassuring feature, since as we noted in the text, many of these models most probably do not 
admit a String/M theory embedding. The fact that the Rushbrooke identity is satisfied suggests that at 
least at large Nc the phenomenological models II are still describing a field theory-like behavior. 

The order of the phase transition is n = [ ^°^^ ] , with \x\ denoting the smallest integer greater or 
equal than x. This permits to accommodate systems with critical exponents /3 > 1/2, the field value 
13 — 1/2 corresponding to interactions containing even powers of "q. Another particular case is /3 = 1, 
which appears in models containing (analytic) interactions if' (i.e. ao = 3). These models undergo 
third order phase transitions. This type of behavior can also be incorporated in the Landau-Ginzburg 
framework in the absence of an 77 -s-)- —77 symmetry. 

We have also studied the conductivity in models obtained by deformations of the HHH model by terms 
0(77*). The conductivity reproduces essentially the same basic behavior as in the HHH model, namely 
a gap a,t Lo < u)g and at high frequencies it approaches the plateau of the normal phase. However, in 
some models the conductivity exhibits resonance peaks in frequency, which are similar to those that have 
previously appeared in the 4+1 dimensional model of [16 . Since the conductivity is directly related to 
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the density of energy eigenstates of the charge carriers, the presence of peaks gives interesting information 
about the energy levels of the condensed matter system. As T — >■ 0, these peaks become very narrow 
and more pronounced. We have explained (following the idea of [2 3) ) how the number of peaks increases 
with the height of the effective potential in a related Schrodinger problem. The examples considered in 
section 5 suggest that there might be models in the family parametrized by the couplings G(ry), U{ri) 
and J{r]) where the potential barrier tends to infinity as T — > 0, thus leading to a hard gap and an 
infinite number of peaks. This would lead to the intriguing possibility that the dual condensed matter 
system would have a density of states given by an infinite sum of delta functions at different frequencies. 
Such behavior is more typical of systems with an external magnetic field, which breaks the Fermi surface 
leading to a discrete structure for the density of energy eigenstates given by the Landau levels. Since in 
the present case there is no magnetic field, in such scenario the discretization of the spectrum might be 
due to some strong coupling effect. Another, possibly more likely, scenario is that, with the back-reaction 
effects incorporated, all models may always have finite number of peaks, that is, any model with finite q 
would have a finite Vmax at T = 0. Therefore, as T — )■ 0, the conductivities would assume a form similar 
to the T = 0.24 curve in fig. 3a, namely no hard gap but extremely suppressed at low temperatures, a 
few peaks and then the DC plateau at sufficiently high frequencies. Even within this more conservative 
scenario, the presence of this finite number of sharp resonances and the gap to the DC plateau becoming 
wider at lower temperatures is quite striking and begs for a condensed matter explanation. Investigating 
this problem in detail requires going to very small temperatures where our numerical approximation with 
no back-reaction is not reliable. 

In the second part of the paper, we investigated the Hall effect induced by a generalized theta term 
of the form 0(r?) F A F. Taking couplings with 0(0) = 0, the effect appears only in the condensed 
phase. The direct and Hall conductivities exhibit quite different features. The direct conductivity has 
a gap, i.e. it is exponentially small at w < Wg, and a delta function at w = representing the infinite 
DC conductivity. The Hall conductivity has no gap and, at w = 0, it has a no n- vanishing, finite value. 
We have also studied the temperature dependence of this finite DC Hall conductivity, which goes to zero 
as T Tc with a critical exponent that depends on 6(77). This type of Hall effect (also studied in a 
different context in [13])) does not require an external magnetic field. Since there are condensed matter 
examples where the Hall effect takes place without any external magnetic field (two examples are the 
topological Hall effect [IS] and the anomalous (spontaneous) AC Hall conductivity ^7\, the latter being 
a manifestation of time-reversal symmetry breaking and emerging in a superconducting state) , it would 
be interesting to see if the generalized "theta" interaction might effectively capture the physics of these 
real systems. 
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